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Material-based two-level systems (TLSs), appearing as defects in low-temperature devices includ¬ 
ing superconducting qubits and photon detectors, are difficult to characterize. In this study we 
apply a uniform dc-electric field across a film to tune the energies of TLSs within. The film is 
embedded in a superconducting resonator such that it forms a circuit quantum electrodynamical 
(cQED) system. The energy of individual TLSs is observed as a function of the known tuning field. 
By studying TLSs for which we can determine the tunneling energy, the actual pz , dipole moments 
projected along the uniform field direction, are individually obtained. A distribution is created with 
60 Pz- We describe the distribution using a model with two dipole moment magnitudes, and a fit 
yields the corresponding values p = pi = 2.8±0.2 Debye and p = p 2 = 8.3±0.4 Debye. For a 
strong-coupled TLS the vacuum-Rabi splitting can be obtained with pz and tunneling energy. This 
allows a measurement of the circuit’s zero-point electric field fluctuations, in a method that does 
not need the electric-field volume. 

PACS numbers: 03.67.-a, 03.67.Pp, 33.15.Kr, 66.35.+a 


Dielectric two-level systems (TLSs) have attracted the 
attention of the quantum computing community ever 
since they were identified as a major source of decoher¬ 
ence in superconducting qubits [1]. Subsequent studies 
found that TLSs were also a performance-limiting source 
of noise in photon detectors used for astronomy [2, 3]. 
This motivation has led to quantum characterization of 
TLSs in the tunneling barrier of superconducting qubits 
[1, 4] and both noise and loss characterization in high 
quality superconducting resonator circuits. Thus TLSs 
are found as defects in various dielectric structures: de¬ 
posited insulating films [1, 5, 6], Josephson Junction (JJ) 
tunneling barriers [1, 7, 8], imperfect interfaces between 
superconducting films with crystalline substrates [9], and 
the native oxides on materials [10]. Recent modeling 
has predicted possible structures and values for the TLS 
dipole moment [11-13]. While these TLSs are generally 
known to be charged atomic configurations that spatially 
tunnel, their microscopic structure and elemental compo¬ 
sition are generally unknown. 

In a qubit made from an anharmonic oscillator the in¬ 
teraction energy is observed as a spectroscopic splitting 
in the qubit state. However, this quantity is not only de¬ 
pendent on circuit element parameters, but also on the 
thickness of the tunneling barrier [1]. The barrier has 
a thickness variation of 1-2 nm [14], such that this lack 
of accuracy is also present in the electric field amplitude 
and the measured TLS dipole moment. In general, indi¬ 
vidual measurements of TLSs within JJs find the tran¬ 
sition dipole moment ptr = where Pz = \Pz\ 

is the absolute value of the dipole moment projected in 
the field direction, and the TLS tunneling energy over 
the total energy Ao/f is generally unknown. As a re¬ 


sult only a lower bound of Pz is determined from those 
measurements (Ao/f < 1). However, in studies of many 
individual TLSs in the barrier of a JJ, the distribution 
of observations is consistent with a single TLS moment 
magnitude p [1, 15], assuming a standard distribution 
of Ao [16, 17]. Here only the average projected moment 
{Pz) ~ p/is determined from the model fit, and has an 
accuracy limitation caused by the aforementioned vari¬ 
able barrier thickness. 

TLSs have been recently studied with an applied field 
which tilts the TLS potential energies. A strain tun¬ 
ing study with individual TLS observations [18] veri¬ 
fied agreement with the standard model which employs 
double-well potentials [16, 17]. Further work revealed 
TLS-TLS interactions [19]. Finally, a swept electric-field 
study measured nonequilibrium microwave loss which 
was explained by Landau-Zener dynamics and a single 
representative dipole moment (magnitude) [20, 21]. 

In our device, TLSs are coupled to a LC resonator and 
tuned with a uniform and known dc electric field 
This tuning distinguishes it from an earlier study of the 
vacuum Rabi splitting and Glauber state made with a 
TLS [22]. TLSs are chosen for analysis only near an 
observed energy minima, which reveals the tunneling en¬ 
ergies Aq. In our experiment, pz = p\ cos6>| is measured 
for each individual TLS and enabled by that of Aq. This 
leaves the relationship to the dipole moment magnitude 
p uncertain only by the angle 0 between the dipole and 
the known dc-field direction. 60 Pz are individually ex¬ 
tracted. The distribution created from these values is fit 
to a model with two dipole moment magnitudes, p = pi 
and p = P 2 , and thus is related to earlier studies using 
two loss tangent sources [23, 24]. While ptr can be de- 
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FIG. 1. (a) Schematic of LC resonator coupled to TLSs within 
the capacitors C1-C4 of the electric field controlled device. 
The voltage Vtias induces a uniform dc-electric field across 
the capacitors, (b) Optical image of the fabricated tunable 
device. Aluminum appears light and the sapphire substrate 
appears black. 


termined in general from dc-tuned measurements, some 
TLSs also have strong-coupling such that the vacuum 
Rabi splitting is observed with the corresponding cou¬ 
pling frequency g = PtrEuMS By combining a mea¬ 
surement of Aq and g we get a complete measure¬ 
ment of the uniform zero-point electric field fluctuations 
Erms of fho device. As will be discussed below, our film 
contains TLSs related to those from JJ barriers, material 
interfaces, contamination, and dielectrics in general, such 
that the accurate moments can generally be extracted by 
this type of device. 

In this experiment we use a superconducting resonator 
with bridge arrangement of four parallel-plate capacitors 
shown in Fig. 1. The capacitors allow a uniform dc elec¬ 
tric bias Ebias to the TLSs within, creating the transition 
energy 

S = y^A2 + (A + 2p,Ebias)2. (1) 

Here, A is the energy difference between the potential 
well minima. The tunneling dipole can alternatively be 
described with a distance between the wells 2p/q for a 
tunneling charge q. The TLS distribution from the stan¬ 
dard model is dn? oc dAodA/Ao [16, 17, 25-27]. Through 
symmetry the bias lead at voltage Fbias is isolated from 
the ac resonance. The capacitors are formed with two 
aluminum films and a do = 125 nm thick layer of silicon 
nitride (SiNx) in between. Since the bridge capacitance 
contains two capacitors in series, the TLSs within expe¬ 
rience the bias field Tobias = i4ias/(2do). The capacitance 
and a quadrupole spiral inductor produce a resonance at 
6.0 GHz. The SiNx dielectric volume within the capac¬ 
itors H = 78 /im^ is micron-scaled (micro-H) to allow 
strong enough coupling for individual TLS observations. 

The dielectric used in this study has a loss tangent of 
tan do = 7.8 x 10“^ as obtained from the Pz data below 
[28]. The silicon nitride is deposited with PECVD[6]. 
Individual TLS have also been observed using films with 
smaller tan do [22], but the film for this study was selected 


to allow many TLSs to be observed in a small range of 
bias voltage Hbias- The microwave input power was gen¬ 
erally set such that the photon occupancy on resonance 
is ^max — 0-4. This photon occupancy allows sufficient 
signal to noise, and in a similar device is found to re¬ 
solve nearly degenerate resonator-TLS energy levels in 
the presence of other weaker coupled TLSs [29]. 

The micro-H device was measured at a temperature 
T < 25 mK in three data sets from separate cool-downs. 
In each data set the bias voltage hbias was varied by 60 
mV across the dielectric corresponding to an electric field 
range of 240 kV/m within it. The bias line was low-pass 
filtered at room temperature and on the 3 Kelvin stage to 
reduce noise induced to the TLS energies. In addition, a 3 
dB microwave attenuator was installed on the 0.7 K plate 
to thermalize the center conductor of the corresponding 
coaxial cable. Finally, a Cu-powder low-pass filter was 
used at the base temperature stage of the refrigerator to 
absorb Cooper-pair breaking photons. 

Figure 2(a) shows the transmission spectroscopy re¬ 
sults from one of the data sets. The wide dark re¬ 
gion generally corresponds to the broad circuit resonance 
co’c/27r, and a large number of fine features are caused by 
near-resonance interactions with individual TLSs. Us¬ 
ing a first analysis which neglects vacuum splitting ef¬ 
fects, we interpret the fine features as the TLS ener¬ 
gies. In agreement with the standard model for energy 
Eq. 1, some TLSs exhibit a minimum in their energy 
(at fiujc — Aq) versus The energy is also sym¬ 

metric about the minimum in the same quantity, as ex¬ 
pected from the referenced equation. Since E^jias is di¬ 
rectly applied, a fit to the data, as shown in Eig. 2(b), 
produces a fairly direct determination of the individual 
projected dipole moments Pz- Each Pz is extracted with 
a precision estimated at approximately 1 percent (see 
below). A = —PzVi)iastdo when the energy is at a mini¬ 
mum {£ = Aq). At this bias, degnerate wells (DW) are 
created in the TLS double well potential (as well as the 
even-amplitude superposition of the single well states). 

Eig. 3 shows the distribution of 60 Pz extracted us¬ 
ing data from three separate cool downs. The result¬ 
ing distribution of Pz can be interpreted to have two 
peaks, and suggests that we compare it to the standard 
amorphous model of TLSs containing two dipole-moment 
magnitudes, pi and p 2 - In accordance with the random- 
orientation defect model, we write the TLS number den¬ 
sity d^n = dpzdAdAo{PiG{pi - pz)/{piAo) + P 20 (P 2 - 
Pz) /(P 2 A 0 )) in terms of the projected dipole moment Pz^ 
where 0 is the Heaviside step function and Pi is the ma¬ 
terial constant for and i = 1,2 [28]. The range of 
dc-bias creates an observable set of DWs with a distribu¬ 
tion 

dnpw _ K PzdVhias , Pl^jpi - Pz) P2Q{P2 - Pz) s 

dPz Wc do Pi P2 

( 2 ) 

Here the experimental ranges are ^Aq = fiK and 6A = 
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FIG. 2. (a) False-color plot of transmission |*S' 2 i| = \Vout/Vin\ vs. frequency / and bias voltage Vbias- Data is taken at T = 25 
mK and fimax — 0.4. Light copper and black correspond to |aS' 2 i| = 0.56 and |aS' 2 i| = 0.40, respectively. TLSs are observed with 
minima in energy (the tunneling energy Aq) nearly degenerate with the circuit photon energy Tiujc- (b) Fits with energy model 
(Eq. 1) to data in panel (a) are shown which find the projected moment pz and tunneling energy Ao of 30 TLSs. 
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FIG. 3. Distribution of pz as determined from fits similar 
to Fig. 2 (b). An initial fit to the model with two dipole 
moment magnitudes p is shown in circle markers and dashed 
lines. The two-sawtooth pattern follows from the technique 
and a standard model of random TLS orientations. The most 
likely ranges for moments magnitudes, pi = 2.8 ± 0.2 Debye 
and p 2 = 8.3 ± 0.4 Debye, are shown by the markers with 
horizontal error bars. 


Pz^ybiasldo^ where tz is the coupling rate between the res¬ 
onator and transmission line. We see from theory that 
for each pi there should be a number of degenerate-well 
TLS in a bin N, where N (x Pz. This feature is a conse¬ 
quence of how larger Pz will be swept further in energy 
for a fixed domain in Vbias- 

An initial fit of Eq. 2 to the data yields pi = 2.9 De¬ 


bye and p 2 = 8.4 Debye, as shown in the figure, where 
circle markers and the dashed line indicate bin heights. 
The fit to the data is consistent with the model if the 
bin to the right of 2.8 Debye has some counts from the 
Pi distribution while the bin to the left of 2.8 Debye has 
lower counts than the fit value due to Poisson statisti¬ 
cal fluctuations. Using Pi and P 2 from the initial fit, 
Monte Carlo analysis was performed using different in¬ 
put values for pi and p 2 - Simulation datasets drew the 
statistically correct number of Pz in simulation distribu¬ 
tions. Fits with optimization for Poisson statistics [30] 
were then used. These fits created simulation-extracted 
Pi and p 2 which generally differ from input values due to 
Poisson sampling error. The difference was used to find 
the probable range of pi and p 2 for a two-moment ma¬ 
terial. Hundreds of simulated datasets were used to find 
the range of moments consistent with the observed data. 
The analysis gives = 2.8T0.2 Debye andp 2 = 8.3T0.4 
Debye for the moments of the model, where the range in¬ 
dicates a single standard deviation. Thus we find our 
data is consistent with a model of two dipole moments, 
and that with the use of this model the standard devia¬ 
tions of dipole moments are <7%. 

To the knowledge of the authors, this experiment 
marks the first study of a single sample from which two 
TLS moment magnitudes were studied. In applications 
with this material we expect that pi is the more impor¬ 
tant moment for material characterization since TLSs in 
this part of the distribution are much more numerous. In 
an alumina barrier of a JJ qubit a TLS dipole moment of 
3 Debye was found if one assumes the barrier thickness of 
2 nm [1] (TLSs detected in tunneling barriers are depen¬ 
dent on an unknown barrier thickness). In ^o-averaged 
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measurements of bulk Si02 with 0H“ correlated TLSs 
the moment magnitude extracted is 3.3 Debye [31]. pi of 
our nitride film is similar and, from SIMS analysis, our 
films are also known to contain some O and H. Although 
our value of p 2 is large, it is within the range of values 
reported in other measurements of silicon nitride [20]. 
We believe that the DW extraction technique creates the 
potential for quantitative analysis on samples, where dif¬ 
ferent dipole distribution models could be tested with 
statistical analysis. 

Furthermore, since our method uses an insulating (cur¬ 
rent blocking) layer for TLS, it can be extended for use 
to characterize various materials which are important to 
quantum-sensitive devices. Any insulating film can be 
tested using this technique, and this includes films of 
alumina, the material commonly used in JJ tunneling 
barriers. In addition, the technique can accommodate 
multiple layers of TLS-laden films, such that it could 
test surface oxides, states at a superconductor/insulator 
boundary, process contamination at the same boundary, 
the thickness dependence of TLSs in films, etc. Finally, 
our technique could test the dependence of p on an in¬ 
dependent variable, such as hydrogen concentration at a 
host material or interface. 

We next analyze one of the most strongly coupled TLS 
to the resonator, where Fig. 4(a) shows a close-up of the 
DW spectra. Here we clearly see the minima energy of 
the TLS with the vacuum Rabi splitting. The Jaynes- 
Cummings eigenenergies for this system are 

E± = -h{iUc + cjtls) ± 9^ + (^/2)^5 (3) 

where cjtls = Sjfi (Eq. I), and S represents the TLS- 
resonator detuning. Figure 4(b) shows the least squares 
Monte Carlo (LSM) fit to the spectra using the theory of 
their lineshapes Ref. [22, 29]. The eigenenergies from the 
fit determined Pz = 6.0 Debye and coupling g/27r = 753 
kHz are shown as white dashed curves. The fit also yields 
the coherence time T 2 = 2/7tls = 313 ns, related to 
the spectral width. Here, 7tls is the TLS relaxation 
rate which limits coherence because the temperature is 
low and strong coupling to the environment, through the 
cavity, provides enhanced relaxation. For this TLS the 
eigenenergy exhibits a visible splitting {g) such that this 
fit yields a more accurate Pz than extracted from the 
previous method, e.g., the method of Fig. 2(b). However, 
the difference in pz caused by even the strongest cQED 
coupling is small; Pz is only different by 1% between the 
two methods. This shows an example of how Pz extracted 
earlier (for Eig. 3) using a weak-coupling assumption are 
only weakly perturbed. Since g/hi = 1.7 is larger than I 
we are in the strong coupling regime. However, the TLS 
relaxation is large in this case, ^/7 tls = 0.74, showing 
that the strong coupling limit does not rely on good TLS 
coherence [32]. 

This analysis also gives a measurement of the RMS 
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FIG. 4. (a) False-color plot of measured |*S' 2 i| vs. frequency 
/ and bias voltage Vbias for one of the most strongly coupled 
TLSs. The resonator-TLS coupling induces avoided crossings 
dependent on the coupling g. Data is taken at T = 24 mK 
and a maximum photon occupancy of Umax — 0.4. Light 
copper and black correspond to |*S' 2 i| = 0.60 and |S' 2 i| = 0.43, 
respectively, (b) Optimum fit to panel (a) data using a cQED- 
based model (see main text). White dashed curves show the 
extracted low-power eigenenergies F+ and E- of the system. 
g is extracted along with pz and Aq, without the electric 
volume. 

fluctuation field Erms without using the geometrically 
estimated electric-field volume. From the optimum fit 
parameters we calculate Erms = T^g^/Pz^o = 25 V/m 
(where Aq/E ^ 1). On the other hand, the design param¬ 
eters allow us to calculate Erms = (2e^eo^E) = 

21 V/m (using = 6.5 [6] and the bridge capacitors’ 
dielectric volume V), which is 16% smaller than the fit 
determined Erms- Since V was not underestimated and 
Pz is very accurate, the difference is believed to be caused 
by the uncertainty in the measured g. 

In conclusion, a resonator was studied which accu¬ 
rately applies a uniform dc-electric field to TLSs in a film. 
The field allows TLS energies to be analyzed as a function 
of a known field and the direct extraction of 60 projected 
moments Pz. An analysis reveals a good fit to a model 
with two representative dipole moment magnitudes for a 
silicon nitride film, pi = 2.8 ± 0.2 and p 2 = 8.3 ± 0.4 De¬ 
bye. The Pi of our nitride film was comparable to that of 
OH“ found in bulk silicon oxide. With additional data 
this method should allow one to search for additional 
dipole magnitudes (TLS types). Since this study uses a 
deposited insulating film containing TLSs, it could be ex¬ 
tended to analyze other materials and structures (surface 
oxides, dielectric-superconductor interface states, fabri¬ 
cation contamination, interdigital capacitors, etc.). This 
should allow the classification of TLSs by dipole moment 
and the optimization of quantum computing or photon 
detector devices. The cQED measurement also allowed a 
complete measurement of the zero-point fluctuation elec¬ 
tric field. This demonstrates how strong-coupling cQED 
can be achieved with TLSs in a practical energy-tunable 
circuit such that many measurements which are avail¬ 
able for superconducting qubits should now be available 
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for TLSs. 
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